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ABSTRACT 


The  long-period  noise  in  the  20  to  40  second  period  range  limits  the  identification 
level  at  which  the  surface-wave,  body-wave  discriminant  can  be  aoplied  at  the  Large 
Aperture  Seismic  Array  (LASA).  Therefore,  an  investigation  was  made  to  determine 
file  sources  and  properties  of  this  noise.  Only  the  long-period  vertical  array  ~t  LASA 
was  considered. 

Both  conventional  and  high -re solution  frequency -wavenumber  spectra  are  pie- 
sented  for  the  noise,  as  well  as  coherence  results.  These  data  show  that  the  noise  con¬ 
sists  of  two  components.  One  component  propagates  across  the  array  as  fundamental - 
mode  Rayleigh  waves  and  is  known  to  be  caused  by  the  action  of  surf  on  coastlines .  The 
other  component  is  nonpropagating  and  evidence  is  presented  which  indicates  it  is 
caused  by  the  elastic  loading  on  the  ground  by  the  earth's  atmosphere.  This  is  estab¬ 
lished  by  correlating  the  power  of  the  nonpropagating  noise  with  the  power  on  the  micro- 
barograph  sensors  at  LASA. 

It  is  also  shown  that  the  signal -to -noise  ratio  gain  obtained  with  maximum - 
likelihood  processing  relative  to  that  obtained  with  beamforming  for  the  long-period 
noise  present  at  LASA,  will  not  be  substantial  unless  it  can  be  shown  that  significant 
amounts  of  propagating  noise  power,  relative  to  total  noise  power,  are  present.  The 
results  at  LASA  indicate  that  such  large  amounts  of  propagating  noise  power  are  rarely 
to  be  observed . 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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I. 


INTRODUCTION 


The  Large  Aperture  Seismic  Array  (LASA)  located  in  eastern  Montana  consists 
of  21  subarrays  of  25  short-period  (SP)  vertical  seismometers  rs  indicated  in  Fig.  1. 

At  the  center  of  each  subarray  there  is  a  three -component  set  of  long-period  (LP) 
seismometers  oriented  in  the  vertical  (Z),  north-south  (NS),  and  east -west  (EW) 
directions.  The  primary  function  oi  LASA  is  to  provide  data  for  facilitating  the  dis¬ 
crimination  between  earthquakes  and  underground  nuclear  detonations. 

It  has  been  found  that  a  useful  discriminant  for  distinguishing  between  natural 
seisrr  *n  events  and  underground  nuclear  explosions  is  based  on  the  relationship  between 
the  surface -wave  magnitude  (M  )  and  the  body -wave  magnitude  (m  ).  1  5  The  surface- 

S  D 

wave  magnitude  M  is  based  on  the  amplitude  of  Rayleigh  waves  with  periods  of  about 
s 

20  seconds  and  is  computed  as^ 


M  =  log  A  —  log  B 
s 

where  A  is  the  ground  amplitude  of  Rayleigh  waves,  in  millimicrons,  with  periods  of 
about  20  seconds  recorded  on  I.PZ  seismometers,  -  log  B  is  a  parameter  which  depends 
on  epicentral  distance  and  is  given  by 


-log  B  =  1.656  log  (A)  -  1. 182 
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wnere  A  is  epicentral  distance  in  degrees.  The  body-wave  magnitude  m^  is  based  on 
the  amplitude  of  SP  waves  recorded  at  teleseismic  distances  and  is  computed  according 
to  the  formula^ 


mfa  =  log  (w/T)  -t 

where  w  is  the  maximum  zero -to -peak  ground  amplitude,  in  millimicrons,  of  the  first 
three  to  four  cycles  of  the  P-vave  recorded  on  vertical  component  SP  seismometers, 

T  is  the  period,  in  seconds,  of  the  observed  short -period  cycle  and  Q  is  a  parameter 
which  depends  on  epicentral  distance  and  focal  depth  and  is  tabulated  in  Fig.  5  of 
Reference  7. 

An  experiment  was  performed  to  determine  the  effectiveness  cf  a  single  LASA 
in  using  the  Mg  -  m^  discriminant  at  teleseismic  distances  and  the  results  have  been 
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given  recently.  Tne  results  of  this  experiment  were  very  encouraging  since  a  perfect 
separation  was  found  to  exist  between  the  Mg  —  m^  characteristic  for  earthquakes  from 
four  tectonic  regions  of  the  earth  and  that  for  presumed  underground  nuclear  exp' op  ions 
from  the  Central  Asian  region.  These  results  indicate  that  if  the  Rayleigh  wave  of  an 
event  can  be  detected  along  with  the  P-wave,  ther  it  is  possible  to  distinguish  whether 
the  source  is  an  earthquake  or  underground  nuclear  explosion  with  very  little,  if  any, 
error. 

Unfortunately,  the  Rayleigh  wave  of  an  event  is  not  as  readily  detectable  as  the 
P-wave  of  the  event.  That  is  to  say,  the  signal -to-noise  ratio  for  the  P-wave  on  the 
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vertical -component  SP  seismometer,  at  about  1  Hz,  is  usualiv  considerably  better  than 
that  for  the  Rayleigh  wave  on  the  vertical -component  LP  aeismometer  at  0.04  Hz.  The 
signal -to -noise  ratio  gain  provided  by  LASA  for  the  P-wave  by  means  of  beamforming 
is  about  20  db.  Similarly,  the  signal -to-noise  ratio  gain  provided  by  LASA  for  the 
Rayleigh  wave  by  means  of  beamforming  is  about  11  db,  and  there  is  an  additional 
8  db  gain  from  matched  filtering,  for  a  total  of  about  19  db.  Thus,  the  signal -to-noise 
ratio  of  the  processed  P-wave  is  still  considerably  better  than  that  of  the  processed 
Rayleigh  wave.  In  other  words,  after  all  of  the  signal  processing  available  at  LASA 
has  been  performed,  it  is  easier  to  detect  the  P-wave  above  the  background  noise  level 
than  it  is  to  detect  the  Rayleigh  wave.  In  this  sense,  it  is  the  background  LP  noise  level 
in  the  20  to  40  second  period  range  which  is  limiting  the  identification  level  at  which  the 
powerful  Mg  -  mfa  discriminant  can  be  applied.  It  is  for  this  reason  that  a  detailed 
mvestigation  of  the  nature  of  the  background  LP  noise  was  made  in  the  20  to  40  second 
period  signal  band. 

It  was  found  that  the  noise  level  on  the  EW  and  NS  LP  seismometi.rs  at  LASA  was 
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usually  about  10  db  higher  than  the  Z  component.  It  is  for  this  reason  that  only  the  Z 
components  were  c  usidered.  This  also  tends  to  simplify  the  interpretation  of  the  data 
and  reduces  the  amount  of  data  processing  required. 

The  amount  of  noise  introduced  by  the  LP  system  was  also  determined.  It  was 
found  that  the  LP  system  contributes  a  negligible  amount  of  noise  power. 
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The  LP  noise  was  analyzed  by  measuring  its  frequency -wavenumber  spectrum 
with  both  conventional  and  high -resolution  frequency -wavenumber  measurement  pro¬ 
grams.  A  detailed  description  of  these  programs  is  given  in  Appendix  A  and  B.  These 
results  indicate  tha<-  the  noise  can  be  considered  as  composed  of  two  components,  a 
component  which  propagates  across  the  LASA  and  a  nonpropagating  component  which 
is  found  to  be  incoherent  over  spatial  lags  greater  than  7.  5  km.  The  propagating  com¬ 
ponent  of  the  LP  noise,  in  the  20  to  40  second  period  range,  is  known  to  be  caused  by 
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the  action  of  surf  on  coastlines.  However,  the  source  of  the  nonpropagating  component 
of  the  noise  is  not  nearly  as  well  established.  In  the  present  v/ork  considerable  evidence 
will  be  presented  to  indicate  that  the  nonpropagating  component  of  the  noise  is  caused 
by  the  elastic  loading  on  the  earth  by  the  atmosphere.  This  result  was  obtained  by 
utilizing  the  microbarograph  sensors  available  at  LASA. 

As  a  consequence  of  the  preceding  results,  the  frequency-wavenumber  structure 
of  the  noise  can  be  determined  quite  accurately.  Using  this  structure  it  is  possible  to 
predict  the  processing  gains  to  be  expected  by  a  multichannel  filtering  method  such  as 
the  maximum -likelihood  method.  These  processing  gains  are  computed  and  presented 
subsequently  and  are  shown  to  agree  with  the  results  obtained  previously  by  processing 
actual  LP  noise  data,  cf.  Reference  8. 
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II. 


NOISE  INTRODUCED  BY  LONG -PERIOD  SYSTEM 


The  first  step  in  the  investigation  of  LP  noise  at  LASA  was  to  determine  Low 
much  noise  was  introduced  by  the  LP  system  shown  in  Fig.  2.  In  this  section  both  the 
EW  and  NS  components  were  considered  as  well  as  the  Z  component.  The  transfei 
function  of  the  LP  system  is  shown  in  Fig.  3.  In  order  to  determine  the  amount  of 
noise  introduced  by  the  LP  system  the  seismometer  mass  was  locked  on  the  EW,  NS, 

Z  components  at  site  D2.  The  spectra  of  these  components,  as  well  as  those  at  sites 
C2  and  Dl,  were  measured  and  the  results  are  shown  in  Fig’.  4.  The  frequency  resolu¬ 
tion  used  in  this  measurement  was  0.013  Hz  and  a  Hanning  window  was  employed.  ^  It 
is  seen  from  Fig.  4  that  the  noise  introduced  by  the  LP  system,  as  given  in  the  spectra 
for  site  D2f  is  about  10  to  30  db  lower  than  that  of  sites  C2,  Dl  in  the  0.025  to  0.05  Hz 
LP  signal  frequency  band.  Since  these  latter  two  sites  can  be  considered  as  having  an 
average  noise  level,  we  may  conclude  that  the  noise  introduced  by  the  LP  system  is 
negligible. 

A  similar  test  was  conducted  to  determine  how  much  of  the  LP  system  noise  is 
due  to  the  LP  system  exclusive  of  the  LP  seismometer.  Towards  this  end  the  EW,  NS, 

Z  seismometers  at  site  D2  were  replaced  with  a  50  K  resistance.  The  spectra  of  all 
three  components  at  site  D2  were  measured  in  the  manner  just  described,  and  are  shown 
in  Fig.  5  along  with  the  results  for  sites  C2  and  Dl.  The  spectrum  for  the  EW  compo¬ 
nent  at  site  C2  is  anomalously  low  due  to  anomalous  signal  on  this  seismometer  at  the 
time.  I*-  is  observed  from  Fig.  5  that  the  noise  levels  a:  site  D2  are  comparable  to 
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those  obtained  during  the  locked  mass  test,  as  seen  in  Fig.  4.  Therefore,  it  may  be 
concluded  that  most  of  the  noise  introduced  by  the  LP  system  is  due  to  the  LP  system 
exclusive  of  the  seismometer. 
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III.  FREQUENCY-WAVENUMBER  STRUCTURE  OF  LONG -PERIOD  NOISE 


The  frequency -wavenumber  spectrum  of  the  LP  seismic  noise  provides  the 
information  about  the  distribution  of  noise  power  with  frequency  and  for  a  fixed  fre¬ 
quency  reveals  the  velocity  and  direction  of  propagation  of  the  noise,  cf.  reference  11 
and  Appendix  A.  Both  conventional  and  high -resolution  frequency -wavenumber  spectra 
of  the  LP  seismic  noise  were  measured.  The  method  of  measurement  is  discussed  in 
detail  in  Appendix  A  and  B.  The  conventional  frequency-wav^.tumber  spectrum  measure¬ 
ment  program  is  due  to  R.  T.  Lacoss.  12  A,  mentioned  previously,  only  the  Z  compo¬ 
nents  of  the  array  will  be  considered.  All  of  the  sensors  at  LASA  were  used  to  measure 
the  frequency -wavenumber  spectrum  except,  for  those  sensors  which  yielded  anomalous 
traces.  The  frequency  resolution  employed  was  0.  01  Hz,  a  Bartlett  window10  was  used, 
and  a  direct  segment,  or  block  averaging,  method  was  employed  with  the  averaging 
performed  over  36  blocks.  This  leads  to  72  degrees  of  freedom  or  90%  confidence  limits 
of  about  ±  i.  2  db.  ^  The  sampling  rate  of  the  data  is  5  Hz  but  a  decimation  factor  of  5 
was  employed  so  that  a  sampling  ra.  a  of  1  Hz  was  used  in  the  measurement.  The  total 
length  of  data  employed  in  the  measurement  was  one  hour. 

The  results  of  both  the  conventional  and  high -resolution  frequency -wavenumber 
programs  are  shown  in  Figs.  6-8  for  three  different  noise  samples  taken  on  26  February 
1967,  7  April  1967  and  26  January  1967,  respectively,  and  for  three  frequencies,  0-03, 
0.04  and  0.05  Hz.  These  figures  show  that  the  conventional  and  high -resolution  results 
are  in  agreement  as  both  methods  tend  to  show  strong  peaks  occurring  at  the  same 
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wavenumber  In  each  progrc  .  .  However,  the  high-resolution  method  delineates  the 

frequency -wavenumber  spectrum  much  more  clearly  than  the  conventional  method.  This 

is  demonstrated  quite  well  in  Fig.  6c,  which  shows  the  wavenumber  structure  of  the 

propagating  noise  to  be  an  arc.  This  is,  of  course,  exactly  what  would  be  expected, 

since  the  dispersion  curve  of  the  propagating  seismic  noise  at  LASA  has  been  measured 
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and  found  to  correspor  1  to  that  of  a  fundamental  mode  Rayleigh  wave.  This  implies 
that  at  a  given  period  the  phase  velocity  of  the  propagating  noise  at  LASA  must  be 
constant,  independent  of  the  sources  of  the  noise,  and  thus  its  wavenumber  structure 
must  consist  of  an  arc,  or  arcs,  whose  extent  corresponds  to  the  range  of  the  azimuths 
of  the  noise  sources. 

An  interesting  example  is  provided  by  the  26  January  1967  noise  sample,  for 
which  the  result3  are  given  in  Fig.  8.  In  particular,  Fig.  8c  shows  a  360  degree 
azimuthal  spread  for  the  wavenumber  structure  with  a  variable  energy  density  along 
this  circle. 

The  preceding  results  show  that  the  LPZ  seismic  noise  consists  of  two  components, 
a  component  which  is  propagating  across  the  array  as  fundamental  -mode  Rayleigh  waves 
and  a  component  which  is  nonpropagating  or  incoherent.  The  nonprojagating  component 
has  been  found  to  be  incoherent  over  distances  greater  than  7.  5  km,  by  measuring  the 
coherence  between  LFZ  sensors  during  time  periods  when  the  propagating  noise  was 
absent  and  noting  that  the  coherence  was  down  to  the  level  corresponding  to  incoherent 
noise.  A  description  of  this  coherence  measurement  program  is  given  in  Appendix  A. 
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A  simulation  was  performed  by  computing  conventional  and  high -resolution  frequ_nc.y 
wavenumber  spectra  for  theoretical  models  of  the  noise  which  correspond  to  f'  ore- 
ceding  description  of  the  noise.  These  simulation  results  were  in  excellent  agreement 
with  those  obtained  for  the  actual  noise  data. 

It  is  quite  important  to  know  the  relative  amounts  of  the  two  kinds  of  noise. 
Towards  this  end,  the  coherence  w~s  measured  for  spatial  lags  between  7.5  and  20  km. 
Since  the  wavelengths  of  the  propagating  noise  are  on  the  order  of  75  to  90  km,  these 
lags  correspond  to  a  small  fraction  of  a  wavelength  and  this  component  of  the  noise  can 
be  considered  to  be  perfectly  coherent  over  these  spatial  lags.  Hence,  any  measured 
loss  of  coherence  over  these  small  spatial  lags  must  be  due  to  the  nonpropagating  com¬ 
ponent  of  the  LPZ  noise.  The  relative  amount  of  power  in  the  propagating  component, 
relative  to  the  total  power  was  obtained  by  averaging  together  the  coherencies  measured 

over  the  small  spatial  lags.  In  addition,  the  bias  introduced  >n  the  measurement  was 

l? 

taken  into  account  by  using  the  tables  of  Amos  and  Koopmans  "  with  their  degrees  of 
freedom  parameter  set  equal  to  the  number  of  blocks  used  in  the  coherence  measure¬ 
ment,  cf.  Fig.  9.  It  is  tacitly  assumed  that  the  noise  power  is  the  same  in  each  LPZ 
sensor.  This  is  not  actually  true,  but  is  a  good  approximation  since  the  noise  power  in 
any  sensor  does  not  usually  vary  from  the  average  noise  power  taken  over  all  of  LASA 
by  about  2-3  db. 

The  relative  amounts  of  noise  power  were  measured  in  this  manner  at  0.03, 
0.04,  and  0.05  Hz  and  the  results  were  then  summed  over  these  frequencies  to  obtain 
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an  approximation  for  the  total  and  relative  amounts  of  noise  power  in  the  0.025  to  0.  05 
Hz  LP  signal  band.  The  results  for  the  total  power  in  this  band  were  checked  independ¬ 
ently  for  two  noise  samples  by  prefiltering  the  data  with  a  sharp  cutoff  400  second  band¬ 
pass  convolutional  filter,  0.025  to  0.05  Hz,  and  then  measuring  the  noise  power  by 
summing  the  squares  of  the  noise  samples.  The  two  methods  gave  results  which  were 
in  agreement  within  1.  5  db.  The  results  of  this  experiment  are  shown  in  Fig.  10  for 
31  noise  samples  extending  in  time  over  a  period  of  about  14  months.  The  results  of 
Fig.  10  show  that  there  is  a  considerable  spread  in  the  total  LPZ  noise  power,  of  about 

12  db.  This  means  that  the  identification  threshold  at  LASA,  using  the  M  -  m, 

s  b 

discriminant,  will  be  variable  depending  on  the  amount  of  background  noise  power  that 
is  present. 

A  cumulative  probability  distribution  was  measured  for  the  ratio  of  nonpropagating 
to  total  noise  power  and  the  results  are  shown  in  Fig.  11.  It  is  seen  from  this  figure 
that  50  percent  of  the  time  this  ratio  is  greater  than  0.  4.  This  fact  has  important  impli¬ 
cations  in  the  consideration  of  array  processing  methods  for  obtaining  signal -to -noise 
ratio  gain,  as  will  be  seen  subsequently. 

A  histogram  o'  noise-source  azimuths  observed  at  LASA  was  measured  at  the  2t> 
second  period  for  the  31  noise  samples  mentioned  previously,  using  the  measured  frequency 
wavenumber  spectra.  This  histogram  is  shown  in  Fig.  12.  It  is  observed  that  much  of  the 
time  the  propagating  noise  tends  to  come  from  the  northwest  and  northeast  directions. 


i 
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IV.  COHERENCE  BETWEEN  LONG -PERIOD  SEEMIC  NOISE  AND 

MICROBAROGRAPH  SENSORS 

It  is  important  to  establish  the  origin  of  the  nonpropagating  LPZ  seismic  noise. 

9 

)*revious  results  due  to  Haubrich  and  MacKenzie  indicate  that  this  noise  might  be  due 
to  the  elastic  loading  on  the  earth  by  the  atmosphere.  Thus,  an  effort  was  made  to  de¬ 
termine  the  coherence  bc<ween  the  LPZ  seismic  noise  and  the  atmospheric  fluctuations 
as  recorded  on  microbarograph  sensors  at  LASA. 

The  LP  seismometers  are  located  in  sealed  metal  tanks  embedded  in  the  floor  of 
an  underground  concrete  structure,  known  as  the  LP  vault.  Each  tank  is  sealed  and  tested 
for  a  leakage  rate  time  constant  of  at  least  eight  hours.  Thus,  any  atmospheric  buoyancy 
effects  on  the  mass  of  the  LP  seismometer,  in  the  20  to  40  second  period  range,  have 
been  effectively  eliminated. 

At  the  time  of  the  present  experiment  there  were  five  microbarograph  installations 
at  LASA,  each  located  near  the  center  of  a  subarray  and  usually  less  than  a  fev.  hundred 
feet  from  the  LP  seisn  ometers.  In  some  cases,  such  as  sites  E3,  B1  and  B4  the  micro¬ 
barograph  was  actually  located  in  the  same  vault  as  the  LP  seismometers.  The  character¬ 
istics  of  the  microbar ographs  are  given  in  Table  I.  The  difference  between  the  two  micro¬ 
barographs  at  site  Ad  is  that  one  of  them  had  various  types  of  wind  filters  while  the  other 
used  a  single  type  of  wind  filter,  namely  a  linear  pipe  array  with  orifices  spaced  10  feet 
apart.  Cor.. parable  results  were  obtained  with  either  of  these  two  microoarographs. 

The  direct  segment,  or  block  averaging,  method  described  in  Appendix  A  was 
used  to  measure  the  coherence  between  LP  seismometers  and  the  microbarographs  and 
also  to  measure  the  spectra  of  these  sensors.  The  frequency  resolution,  frequency  window, 
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CHARACTERISTICS  CF  MICROBAROGRAPHS 
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number  of  blocks,  anti  decimation  factor  are  the  same  as  that  used  previously  in 
Section  111  in  the  measurement  of  frequency -wavenumber  spectra. 

Typical  spectra  for  the  LPZ  sensor  at  site  A0  and  one  of  the  microbarographs  at 
site  A0  are  shown  in  Fig.  13.  The  dates  and  times  of  the  noise  samples  used  to  measure 
the  coherence  are  listed  in  Table  II.  It  should  be  noted  that  only  9  of  the  17  noise  samples 
showed  any  coherence  between  an  LFZ  seismometer  at  any  site  and  its  corresponding 
microbarograph.  It  is  also  worth  noting  that  the  noise  samples  were  picked  so  as  to  have 
a  reasonably  large  amount  of  nonpropagaring  noise  in  the  20  to  40  second  period  range. 

Two  examples  in  which  coherence  was  measured  at  site  A0  are  given  in  Figs.  14a  and  b, 
where  the  latter  figure  is  more  typical  of  the  behavior  of  the  coherence  when  any  coherence 
is  measured  at  all.  In  Fig.  14a  the  coherence  at  0.03  Hz  is  quite  high,  about  0.6  and 
drops  to  about  0.5  and  0.4  at  0.04  and  0.05  Hz,  respectively.  The  ratio  of  nonpropagating 
noise  to  the  total  noise,  for  this  noise  sample,  was  measured  as  0.  9  at  0.03,  0.04  and 
0.05  Hz.  The  measured  coherence  is  not  as  high  as  this  at  these  frequencies,  but  is  the 
largest  that  was  ever  measured  in  all  of  the  17  noise  samples.  The  results  of  Fig.  14b 
are  more  typical  in  that  whenever  coherence  is  measured,  it  tends  to  be  Irgh  at  0.02  Hz 
and  then  drops  to  the  level  for  incoherent  noise  at  0.03  to  0.05  Hz.  The  95  percent  confi¬ 
dence  limits  for  this  level  for  incoherent  noise  can  be  obtained  from  Fig.  9  a r  0. 02  to 
0.  3.  The  ratio  of  nonpropagating  to  total  noise  for  the  noise  sample  used  in  Fig.  14b  was 
measured  as  0.60,  0.45,  and  0.35  at  0.03  ,  0.04  and  0.05  Hz,  respectively.  It  is  seen 
that  the  coherence  never  reaches  a  level  compatible  with  the  amount  of  non  propagating 
noise.  Thus,  we  may  conclude  that  there  is  relatively  little  co!  erence  between  the  LFZ 
noise  and  the  microbarograph  signals  in  the  20  to  40  second  period  range. 


13 


TABLE  II 


NOISE  SAMPLES  USED  TO  DETERMINE  COHERENCE  BETWEEN 
LPZ  SEISMOMETERS  AND  MICROBAROGRAPHS 


DATE 

START 

TIME 

(GMT) 

28  May  67* 

03:40:00 

3  June  67* 

10:47:00 

1  July  67* 

19:50:00 

10  July  67* 

02:30:00 

15  July  67* 

02:50:00 

29  July  67* 

07:52:00 

23  Aug  67* 

17:20:00 

2  Sept  67 

01:43:00 

30  Sept  67* 

01:12:00 

1  Oct  67 

00:33:00 

1  Nov  67 

13:53:00 

24  Nov  67 

16:25:00 

31  Dec  67* 

04:05:00 

5  Jan  68* 

15:30:00 

5  Jan  68 

23:30:00 

6  Jan  68 

19:30:00 

7  Jan  68 

04:45:00 

*  Denotes  that  coherence  was  observed  between  LPZ  seismometer  and  microbarograph. 

EACH  NOISE  SAMPLE  IS  ONE  HOUR  LONG. 
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An  attempt  was  made  to  correlate  the  amount  of  nonpropagating  LPZ  seismic 
noise  power  and  the  amount  of  microbarograph  noise  power  in  the  20  to  40  second 
period  range.  These  quantities  were  measured  at  site  AfJ  for  the  17  noise  samples 
given  in  Table  II  and  the  results  are  shown  in  Fig.  15.  This  figure  shows  that  there  is 
a  definite  trend  for  the  power  level  on  the  LPZ  seismometer  to  increase  when  the  power 
level  on  the  microbarograph  increases.  Thus,  this  is  evidence  that  the  nonpropagating 
LPZ  noise  is  caused  by  the  elastic  loading  on  the  earth  by  the  atmosphere.  The  non- 
propagating  seismic  noise  may  be  caused  by  the  cumulative  effects  of  many  independent 
atmospheric  pressure  fluctuation  cells  within  a  radius  of  a  few  km  of  the  LPZ  sensor 
and  this  may  be  the  reason  for  the  low  coherence  between  this  noise  waveform  and  a 
single  microbarograph  record  obtained  at  the  same  location. 
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V. 


COMPUTATION  OF  ARRAY  PROCESSING  GAIN  FOR  THEORETICAL  MODEL 
OF  SIGNAL  AND  NOISE 


As  a  consequence  of  our  preceding  results  it  is  possible  to  give  a  rather  good 
representation  for  ':he  structure  of  the  LPZ  seismic  noise,  in  the  20  to  40  second  period 
range,  in  frequency -wavenumber  space.  It  has  been  established  chat  the  noire  consists 
of  two  components.  One  component  is  propagating  across  the  LASA  as  a  fundamental- 
n  ode  Rayleigh  wave  and  its  wavenumber  structure,  at  a  particular  frequency,  consists 
of  an  arc,  or  possibly,  arcs.  The  radius  of  the  arc,  ,  determines  the  phase  velocity 
of  the  noise,  v,  by  means  of  the  equation  v  =  f/k,  where  f  is  the  frequency.  The  angular 
extent  of  the  arc,  or  arcs,  is  determined  by  the  azimuthal  distribut*  in  of  the  noise. 

The  other  component  of  the  noise  is  noupropagating  and  incoherent  and  its  coherence 
structure  consists  of  a  Kronecker  delta  function,  i.e. ,  the  coherence  is  seme  value, 
between  zero  and  unity,  for  zero  spatial  lag  and  is  zero  otherwise. 

It  is  important  to  predict  the  signal-to-noise  ratio  gain  achievable  by  various 
forms  of  array  processing.  If  the  wavenumber  structure  of  the  noise  is  known,  at  a 
given  frequency,  then  for  a  fixed  seismic  array  geometry  it  is  possible  to  compute  the 
signal-to-noise  ratio  gain  of  a  particular  form  of  array  processing.  This  computation 
was  performed  using  the  array  at  LASA  formed  from  the  A,  C,  D,  E  and  F  rings,  and 
the  details  of  the  computation  are  presented  in  Appendix  C.  The  array  processing 
methods  considered  were  the  maximum -likelihood  and  del  ay -and -sum,  or  beamforming, 
methods,  cf.  reference  14.  The  results  of  the  computation  are  shewn  in  Fig.  16  which 
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plots  the  gain  of  maximum -likelihood  relative  to  beamforir.irg  vs  azimuth  of  the  noise 
source  and  the  extent  of  the  arc.  The  wavenumber  structure  of  the  noise  is  assumed 
to  consist  of  a  single  arc  and  the  azimuth  of  the  noise  source  is  defined  as  the  azimuth 
to  the  center  of  this  arc.  The  array  is  assumed  to  be  steered  for  an  event  from  due 
north  whose  phase  velocity  is  3.7  km/sec,  the  frequency  is  0.04  Hz  and  the  phase 
velocity  of  the  noise  is  also  3.7  km/sec.  Results  for  other  event  steering  parameters 
yielded  comparable  results.  The  ratio  of  incoherent  to  total  noise  power,  R,  is  0.1, 

0.  3  and  0. 5  in  Figs.  16a,  b,  and  c,  respectively.  These  results  show  that  when  R  =  C.  1 
there  is  significant  gain  of  maximum -likelihood  relative  to  beamforming,  for  a  wide 
range  of  noise  source  azimuths  and  arc  extents.  However,  the  condition  R  =  0. 1  is 
rarely  met  at  LASA  as  can  be  seen  Lrom  Fig.  11.  The  results  in  Figs.  16b  and  c  show 
that  there  is  not  much  gain  for  maximum -likeiihood  relative  to  beamforming.  These 
conditions, when  R  =  0.  3,  0.  5,are  more  typical  of  the  situation  encountered  at  LASA. 

It  should  be  mentioned  that  these  theoretical  results  were  checked  with  results  obtained 
with  actual  data  with  reasonable  agreement  to  within  2  db.  The  B-ring  was  omitted  from 
the  computations  since  the  effect  of  including  it  would  have  been  to  lower  the  gain  of 
beamforming  considerably  and  increase  the  gain  of  maximum-likelihood  only  slightly, 
due  to  the  small  diameter  of  the  B-ring.  When  the  B-ring  is  omitted,  the  amplitude 
gain  of  beamforming  is  approximately  equal  to  ,/NS,  where  NS  is  the  number  of  sensors 
in  the  array,  provided  the  azimuth  of  the  propagating  noise  is  not  within  a  haif  beamwidth 
of  the  azimuth  of  the  event.  Note  that.  NS  =  21  —  4  =  17,  ,/NS  =  5.4,  corresponds 
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to  12.  3  db  of  gain.  This  aids  in  providing  ar.  approximate  absolute  scale  ior  the  gain  of 
maximum  -likelihood  and  beamforrning  methods  when  using  the  results  of  F;g.  16. 

Thus,  unless  it  can  be  established  that  significant  amounts  of  noise  are 
propagating  s  a  fundamental -mode  Rayleigh  wave,  the  form  of  arra^  processing  to 
use  is  beamforming,  because  of  its  simplicity  of  implementation.  The  rnaximum- 
like'ihood  method  is  much  more  complicated  and  does  not  provide  significant  gain  over 
beamforming  unless  R  <  0.  3,  a  condition  which  is  rarely  met  at  LASA.  However,  if 
it  were  possible  to  reduce  the  amount  of  nonpropagating  noise,  then  the  maximum  - 
likelihood  method  would  be  very  useful  in  obtaining  large  signal -to-noise  ratio  gains. 
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VI.  CONCLUSIONS 


The  LPZ  background  seismic  noise  is  limiting  the  identification  threshold  at 
which  the  Mg  —  m^  discriminant  can  be  applied.  This  noise  h.  i  been  investigated  and 
found  to  consist  of  two  components.  One  component  propagates  across  LASA  as  a 
fundamental -mode  Rayleigh  wave,  and  is  caused  by  the  action  of  surf  on  coastlines. 

The  other  component  is  nonpropagating,  or  incoherent  noise,  and  data  have  been  ob¬ 
tained  which  indicate  that  it  is  caused  by  the  elastic  loading  on  the  ground  by  the  earth’s 
atmosphere.  This  component  always  provides  a  significant  contribution  to  the  total 
noise  power.  As  a  consequence,  it  has  been  shown  that  sophisticated  array  processing 
methods  such  as  maximum -likelihood  cannot  provide  significant  gain  relative  to  beam¬ 
forming. 

If  the  nonpropagating  noise  is  caused  by  the  action  of  the  atmosphere  on  the 
ground,  then  the  amplitude  of  this  noise  will  be  determined  by  the  Lame  parameters 
of  the  medium,  cf.  reference  9,  p.  1440.  The  LASA  is  located  on  sedimentary  layers 
in  eastern  Montana,  so  that  the  Lame  parameters  are  relatively  small.  It  is,  therefore, 
desirable  to  construct  any  future  array  on  granite  whose  Lame  parameters  will  be 
larger  than  those  for  sediments.  This  will  lead  to  a  reduction  of  the  power  of  the  non¬ 
propagating  component  of  noise.  In  this  case  there  may  be  an  important  role  to  be 
pla}  xl  by  sophisticated  array  processing  methods  such  as  maximum -likelihood. 
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APPENDIX  A 


DESCRIPTION  OF  CONVEN HONAL  FREQUENCY- WAVENUMBER 
SPECTRUM  MEASUREMENT  PROGRAM 

We  assume  that  the  noise  {N  },  j  =  1,. . .  ,K,  K  is  the  number  of  sensors,  is 
a  wide-sense  stationary  discrete -time  parameter  random  process.  The  covariance 
matrix  of  the  noise  is  given  by 


p  (m-n)  =  E  {N.  N  },  j,k  = 
jk  jm  km  J  ’  ’ 


(1) 


where  E  denotes  expectation.  The  cross -power  spectral  density  is 


fjkM  =  X  pjk(m)  j.k=1«---.K 

J  m=-®  J 


(2) 


and 


TT 


ym>  -  J-  e 

J  -  TT  J 


-imX  dX 

2tt 


(3) 


where  X  =  2tt{T,  f  is  frequency  in  Hz  and  T  is  the  sampling  period  of  the  data  in  seconds. 

If  the  noise  field  is  space  stationary,  then,  for  fixed  X,  f  (X)  depends  only  on  the 

j-o 

vector  difference  x  .  —  x  , ,  where  x  .  is  the  vector  position  of  the  sensor.  It  is  con- 
-j  -l  -j 

venient  to  introduce  a  cross -power  spectral  density  f(X,r)  as 
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f(\,r)  =  f^(X),  if  —  x ^  =  v 

and  is  also  defined  for  all  possible  vector  lags  r  in  a  similar  manner.  That  is,  we  can 
imagine  that  sensors  have  been  placed  at  any  desired  positions  so  that  we  can  define 
f(\,  r)  as  given  aoove.  It  is  now  possible  to  define  the  frequency -wavenumber  spectrum 
of  the  noise  as 


00  CO 


P(\,  k)  =  J  f  f(\,  r)  €l—  —  drxdry 


(4) 


where  k  is  die  vector  wavenumber  and  r  ,  r  are  the  x,  y  components,  respectively, 

x  y 

of  the  vector  r.  We  may  now  use  Eq.  (4)  to  obtain  f  (X)  as 


00  eo 

=  !  J  P(X,k)  €'l-'(-j  dkxdky  j,l=l . K 


(5) 


where  k  ,  k  are  the  x,  y  components,  respectively,  of  the  vector  k,  in  radians /km. 
x  y 

If  the  noise  consists  of  a  unity  amplitude  monochromatic  plane  wave  propagating 
with  a  velocity  v  ,  of  the  form  ^l^2TTloj^  — o  E.)^  j  =  o,  ±  1,  ±  ?.,... ,  where  fQ  is  the 
frequency, kQ  =  2n  fQ  a^,  aQ  is  a  slowness  vector  which  points  in  the  direction  of  propa¬ 
gation  of  the  wave  and  [  |  =  1/v^,  then 


f(X,r)  = 


and 
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dr  dr 
x  y 


P(X,k) 


00  00 


-  J- J 


€i(k-k0) 


=  Mk-k0) 


which  is  a  delta  function  located  at  the  wavenumber  kQ.  It  should  now  be  apparent  how 
P(X,  k)  provides  the  information  concerning  the  spe^d  and  azimuth  of  propagating  seismic 
waves. 

The  program  which  provides  an  estimate  for  P(X,  k)  employs  the  direct  segment, 
or  block  averting,  method  for  reasons  of  computational  efficiency,  as  described  in 
reference  14.  The  number  of  data  points,  L,  in  each  channel  is  divided  into  M  non¬ 
overlapping  blocks  of  N  data  points.  The  Fourier  transform  of  the  noise  data  in  the 
nth  segment,  jth  channel,  and  normalized  frequency  X,  is 


S, 

jn 


l 


m=l 


n  eim>v 

j,  m  +  (n-  1)N 


j  =  1 . K 

n  -  1,..., M . 


(6) 


As  an  estimate  for  f.,  (X)  we  take 
Jk 


M 

-  m  i  yi>  ■ 

J  n=  1 


j, k  -  1 , . .  •  ,  K  . 


(7) 


At  this  point  the  program  performs  a  normalization  by  dividing  f  (X)  by  [f..(X)f  (X)] 

jk  JJ  kJc 

in  order  to  remove  the  effects  of  improper  sensor  equalization.  We  can,  without  any 
loss  of  generality,  ignore  this  step  in  the  ensuing  analysis. 


1/2 
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As  an  estimate  for  P(\,k)  we  take 


P(M)  =  i  S 


It  was  shown  in  reference  14  that  {  f  (X)  }  is  a  nonnegative -definite  matrix  so  that 

A  /S 

P(X,k)  will  be  real  and  nonnegarive.  Ue^’ng  the  result  in  reference  14  for  E  (f  (X))  we 


TT  00  co 


E{P(X,ko)}  -  J  J*  J  P(\,k)  |B(k-ko)r'-  |WN(x->.)|2 

—  TT  —00  —00  ^ 


2  10 
where  |W^.(x)|  is  the  Bartlett  window 


|W  fxll2  =  -  sin  <N/2>  x 

1  NWI  N  sin  (1/2)  x 


and  |B(k)j  is  the  beamforming  array  response  pattern 


boo  =  i  1  . 

pi 


Thus,  E{P(\,k  )}  is  obtained  by  means  of  a  frequency -wavenumber  window 
2  2 

|W^(x-X)|  |  B(k  -ko)|  .  Hence,  P  will  be  an  asymptotically  unbiased  estimate  for 

2 

cP  if  |  WN  (x  —  \)  •  B  (k  —  k  q)|  approaches  a  delta  function  in  such  a  way  that 


J  J  J  IWN(x-X)-  B(k-ko)|2  |  dkx  dk  =  c 


— TT  —oo  —oo 


23 


where  e  is  some  positive  constant.  The  fceamforming  array  response  pattern  for  the 

K  ky 

LPZ  array  at  LASA  is  given  in  Fig.  17  which  shows  -  20  log  |B(k)|  vs  ^  ,  in 
cycles/km. 

A 

Using  the  results  of  reference  14  we  can  compute  the  variance  of  P  as,  assuming 
{N  }  is  a  multidimensional  Gaussian  process, 


vAR[p<x,ko)]  =  {e  ip<x,ko)i>2  +  i  j  J  /  jw,sa%yB4ty 


'-TT— eo  — * 


lWN<‘-x>|2  K  dk*  dkv 


Thus 


VAR  (P(X,ko)]  ^  ±  {E[P(X,ko)]}2,  (  kQ|  t  0 

=  ^  {E[P(X,ko)]}2,  |kj  =  0.  (12) 

A 

The  mean  and  variance  of  P  are  given  by  Eqs.  (9)  and  (12),  respectively. 

We  follow  Blackman  and  Tukey,1^  and  assume  that  P(^,kQ)  is  a  multiple  of  a 
chi-square  variable  so  that  to  establish  confidence  intervals  the  chi-square  distribution 
can  be  used  with  the  number  of  degrees  of  freedom  k  given  by 

k  =  2{E  (P(X,ko))}2/VAR(P(X,ko)] 

=  2M,  jk  |  t  0 

—  o 

=  M,  Ik  I  =  0  . 

'-o' 
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If  M  =  36,  k  =  72  and  the  90  percent  confidence  limits  are  approximately  ±1.2  db,  if 
|k  |  ^  0.  When  |k^|  =  0,  these  limits  are  approximately  ±  1.6db. 

If  the  noise  consists  of  a  plane  wave  propagating  across  the  array,  plus 
incoherent  noise,  then,  if  M  and  N  are  large,  the  spectral  matrix  is  given  by 

fK(X)  =  V(R)  . K*  (i3> 

where 

5^(R)  =  1  ,  )  =  l 

=  1-R,  \t>  , 

R  is  the  ratio  of  incoherent  to  total  noise  power,  kQ  =  2rrfa,  f  is  the  frequency  of  the 
noise,  a  is  the  slowness  vector  which  points  in  the  direction  of  propagation  and  has  mag¬ 
nitude  | at |  =  1/v  ,  v  is  the  phase  velocity  of  the  propagating  component  of  the  noise. 
Hence,  we  have  using  Eq.  (8) 

PO^k  )  =  1  -  R  +  § 

— o  K 

The  frequency -wavenumber  spectrum  measurement  program  displays  contours,  at  a 

A  A  A 

fixed  frequency,  of  - 10  log  [  P(\,k)/Pj^^(\,k)]  vs  V  ^/2n,  k  /2tt,  where  P  v  is  ^ 

A 

maximum  value  of  P.  The  wavenumber  coordinates  are  in  cycles/km.  In  the  present 
case  the  contours  will  tend  to  peak  up  at  the  wavenumber  k  =  k^,  and  if  R  is  small  will 
tend  to  have  contours  similar  to  those  of  the  beam  pattern  B(k),  cf.  Eq.  (11). 


We  close  this  section  by  defining  the  sample  coherence  between  the  and  k1*1 


sensors,  as  used  previously  in  the  measurements  on  noise  data, 


C  (X)  =  - £ - 


1/2 


and  the  coherence  is,  of  course, 


(X)  = 
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APPENDIX  B 


DESCRIPTION  OF  HIGH-RESOLUTION  FREQUENCY -WAVENUMBER 
SPECTRUM  MEASUREMENT  PROGRAM 


The  high-resolution  estimate  for  P(\,k)  is  defined  as 


r  k  . 

’■<>'.«  =  E  q,. 


-i 


(14) 


where  {q  (X)S-  is  the  inverse  of  the  spectral  matrix  {f .  .(X)}.  The  motivation  for  this 
jv  j  K> 


procedure  can  be  given  by  writing  Eq,  (14)  as 
K 


.  M  ,  K 

1  V 


?’(X,k)  =  Yi  a  (X,k)A  (X,k)£  (X)*-'*-!-^  =  M  ?  Z  Aj<X**>  Sjn<X) 
j.t-l  J  1  J 


where 


K 

Z 

A.(X,k)  =  - 

Z  qu(^JE) 

J 


(15) 


and  { q  (X , k)}  is  the  inverse  of  the  matrix  {€*-  ^  f..(X)}.  Thus,  P’(X,k  )  is  the 

j J'v  ^ 

output  of  a  maximum -likelihood  filter,  whose  design  is  determined  by  the  noise  data  and 

is  different  for  each  wavenumber  k0,  which  passes  undistorted  any  monochromatic  plane 

wave  travelling  at  a  velocity  corresponding  to  the  wavenumber  k  and  suppresses  in  an 

o 

optimum  manner  the  power  of  those  noise  waves  travelling  at  velocities  corresponding 
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to  wavenumbers  other  than  k  ,  cf.  reference  l*1  It  should  be  noted  that  the  amount  of 

— o 

computation  required  to  obtain  P'  is  almost  the  same  as  that  to  get  P,  since  only  an 
additional  Hermitian  matrix  inversion  is  required. 

We  now  wish  to  compute  the  mean  and  variance  of  P1.  In  order  to  do  this  it  will 
be  assumed  that  the  weights  A^(X,k)  defined  in  Eq.  (15)  are  not  random.  This  is  a 
simplifying  assumption,  which  is  not  actually  valid,  since  these  weights  are  designed 
from  the  noise  data.  However,  it  does  appear  to  be  a  reasonable  approximation  and 
the  results  obtained  would  be  valid  if  these  weights  were  designed  somehow  independent 
of  the  noise  data.  Using  this  assumption  we  have 


2  dx 
2tt 


dk  dk 
x  y 


(10) 


where 


B'(X,k,  kQ) 


K 

l 

j=i 


A.(X,k„> 


€i(k-k0)*Xj 


(17) 


Thus,  E[P‘(X,ko)]  is  obtained  by  means  of  a  frequency -wavenumber  window 

2 

|  W^(x  -  X)  •  B'  (X,  k,  kQ)  |  .  Hence,  P'  will  be  an  asymptotically  unbiased  estimate  for 

i  i2 

cP  if  j  (x  -  X)  •  B'  (X,  k  ,  k^)  |  approaches  a  three-dimensional  delta  function  in  such 
a  way  that 

iff  |WN<*-X).  B-(A,k,ko)!2  f  dkx  dk  -  c 

— T7  —<»  — oo  J 
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where  c  is  some  positive  number. 


The  variance  of  P'  is,  assuming  {N  }  is  a  multidimensional  Gaussian  process, 


V‘R[P'<\,k„S]  =g{E|P'(X,ko)]}2+i 


—  0 


TT 


J  J  J  P(X,k)B,*(\,ktko)B,(\tk,-k)) 


-jj  -oo  -co 


|  W  (x  —X) | “  p-  dk  dk 
N  1  2tt  x  y 


Thus , 


VAR  [P’(X,k  )]  s  ~  {E  lP'(X,k.  )]  }2,  |k|  4  0 

■  R  {HiP(X,k&)|}2,  |ko|-0. 

A 

The  confidence  limits  for  P'  can  be  obtained  in  a  manner  similar  to  that  for  P  described 
in  Appendix  A. 

It  is  interesting  to  compute  P'  for  the  spectral  matrix  given  in  Eq.  (13).  If  we 
denote  this  matrix  by  F,  then 

F  =  (1-R)  [q'  q  +  1  1 

where  I  is  the  K  x  K  identity  matrix,  q  is  a  l  x  K  row  matrix  whose  1  element  is 
G*-°  -j  and  q'  is  a  K  x  1  column  m  itrix  whose  j  1^  element  is  G  —  j.  Now  we 

have 
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l(l-R)(q,q  + 


t?r  1)1  1  ■  5“- 


li 


k  + 


R 

1  -R 


so  that 


P*(X,k)  = 


R  rR  +  k (1  -  R)] 


2  2R  a 

K  [1-  R+~  -P(X,k)] 


Now  when  k  =  k  ,  P(X,k  '  =  1-R  +  -  and 
-  — o  —o'  K 


P’(X,k  )  =  1-R+  | 
— o  K 


=  P(X,k_o) 


The  high  resolution  frequency -wavenumber  spectrum  measurement  program  displays 

contours,  at  a  fixed  frequency,  of  -  10  log  [P'(X , y/P^j^X.k)]  vs  —  ,  ^  ,  where 

P! -AV  is  the  maximum  value  of  P'.  The  wavenumber  coordinates  are  in  cyclee/km. 
MAX 

It  will  now  be  shown  that  the  wavenumber  resolution  using  )?'  is  higher  than  that 

A 

of  P.  We  assume  the  noise  is  such  that  the  spectral  matrix  is  given  by  Eq.  (13).  In  the 
vicinity  of  k  =  kQ  we  consider  the  contour  for  which  P*(X,k)  =  1  -R,  which  is  still  very 

p 

close  to  the  peak  value  of  1  -R  +  —  ,  since  R  is  smail,  between  zero  and  unity,  and  K 

N. 

is  large,  usually  about  20.  For  these  values  of  k  we  have 


P' (X r k )  -  J  (i~R+  |) 
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so  that  P’  is  already  three  db  down  from  its  peak  value  of  1  —  R  +  —  .  Hence,  the 

K. 

wavenumber  resolution  of  P'  will  be  much  higher  than  that  of  P. 

We  close  this  section  by  noting  that  the  high- resolution  program  described  is 

15  16 

similar  to  c  te  proposed  by  Burg  and  used  by  McCowan  a-td  Lintz.  Although  the 
two  methods  are  different,  it  may  be  shown  that  both  methods  yield  similar  results. 
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APPENDIX  C 


DESCRIPTION  OF  COMPUTATIONS  FOR  OBTAINING  GAIN  OF 
BEAMFORMING  AND  MAXIMUM -LIKELIHOOD  PROCESSING 


We  now  wish  to  present  a  description,  of  the  computations  used  to  obtain  the 
results  presented  in  Seciion  V  concerning  the  gain  of  beamforming  and  maximum  - 
likelihood  processing.  The  noise  is  considered  to  consist  of  two  components,  a  com¬ 
ponent  which  is  propagating  across  the  array  as  fundamental -mode  Rayleigh  waves  and 
a  nonpropagating  component,  as  described  previously.  The  spectral  matrix  of  the  noise 
is  computed  by  means  of  the  following  numerical  integration  over  the  wavenumber 
structure  of  the  noise 


100 


fjk(f)  =  ait(R)€i2TTfi-^rr-k)^_  £  €i2TTfSn-ik 


jk 


100 


n=l 


where  6  (R)  was  defined  previously  in  Appendix  A,  following  Eq.  (13),  f  is  the  frequency 
in  Hz,  B  is  a  slowness  vector  which  points  in  the  direction  of  propagation  of  the  event, 
for  which  the  array  is  being  steered  by  insertion  of  appropriate  time  delays,  whose  mag¬ 
nitude  is  1/V„,  V  is  the  Dhase  velocity  of  the  event  in  km/sec,  r  .  is  the  vector  position 
of  the  jth  sensor  in  km,  r  ^  is  the  vector  position  of  the  ktfl  relative  to  the  jth  sensor  in 
km,  an  is  a  slowness  vector  which  points  outward  from  the  oiigin  in  wavenumber  space 
in  the  direction  of  the  azimuth  9  -  |  +  n-^  and  its  magnitude  is  1/VN,  VN  is  the 

phase  velocity  of  the  noise  in  km /sec,  <b  is  the  angular  extent  of  the  arc  in  wavenumber 
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space  and  B^  is  the  azimuth  to  the  center  of  this  arc.  The  beamforming  gain,  in  db, 
is 


GB  -  -10  1og10(K-2  V  f  <01 

1,  k=  1 


The  gain  of  maximum -likelihood  processing,  in  db,  is 


K 

GML  =  10  log  [  2]  q*  (0  ] 
1U  j,k=l  J 


where  the  matrix  {q*  (f)}  is  the  inverse  of  { f  (f)} .  The  gain  of  ma 
jk  jk 


mum-likel  ihood 


processing  relative  to  beamforming  is  (GML  -  GB). 
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Fig.  1.  General  arrangement  of  the  large  aperture  seismic  array. 
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3.  Long-period  system  transfer  function  (25-second  filter) 
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POWER  DENSITY  (aB  relative  to  tmnVHz  ot  0.04  Hz  ) 
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16.41  00  TO  17:17.20 


Fig.  4.  Power  spectral  density  results  for  locked  mass  test. 
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SEISMOMETERS  AT  SITE  02  REPLACED  BY  50K  RESISTANCE 


Fig.  5.  Power  spectral  density  results  when  seismometers  at  site  D2  are  replaced 
by  50  K  resistance. 
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RESOluT 


olution  frequency-wavenumber 


Continued. 


Conventional  and  high-resolution  frequency-wavenumber  spectra  for 
967  noise  sample. 


Fig.  7.  Continued. 
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Fig.  8.  Conventional  and  high-resolution  frequency-wavenumber  spectra  for 
26  January  1%7  noise  sample. 
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Fig.  8.  Continued. 


SAMPLE  COHERENCE 


Fig.  9.  The  mean  and  95  percent  confidence  limits  for  sample  coherence  using 
direct  segment  method  with  36  blocks. 
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Fig.  10.  Propagating,  nonpropagating  and  total  noise  power  in  the  20  to  40  second 
period  range  vs  time  for  31  noise  samples. 
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Fig.  11.  Measured  cumulative  probability  distribution  function  for  ratio  of  nonpropagating 
to  total  noise  power  in  the  20  to  40  second  period  range  for  31  noise  samples. 
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FREQUENCY  =  0  04  Hz 


Fig.  12.  Histogram  of  noise-source  azimuths  observed  at  the  large  aperture  seismic 
array  at  25  second  period  for  31  noise  samples. 
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Fig.  13.  Power  spectral  densities  for  long-period  seismic  noise  and  microbarograph 
record  at  site  A0. 
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Fig.  14.  Coherence  between  long-period  vertical  seismometer  and  micrebarograph, 
at  site  A0,  for  two  noise  samples. 
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Fig.  15.  Noupropagating  seismic  noise  power  vs^  microbarograph  noise  power  in  the 
20  to  40  second  period  range,  at  site  A0. 
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Fig.  16.  Gain  of  maximum- likelihood  processing  relative  to  beamforming. 
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Fig.  16.  Continued. 
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Fig.  16.  Continued. 
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Fig.  17.  The  beamforming  array  response  pattern  for  the  large  aperture 
seismic  array. 
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